Locality for quantum systems on graphs depends on the number field 
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Adapting a definition of Aaronson and Ambainis [Theory Comput. 1 (2005), 47-79], we call a 
, quantum dynamics on a digraph saturated Z-local if the nonzero transition amplitudes specifying 

the unitary evolution are in exact correspondence with the directed edges (including loops) of the 
digraph. This idea appears recurrently in a variety of contexts including angular momentum, quan- 
£SJ , turn chaos, and combinatorial matrix theory. Complete characterization of the digraph properties 

that allow such a process to exist is a long-standing open question that can also be formulated in 
f"S |' terms of minimum rank problems. We prove that saturated Z-local dynamics involving complex 

| amplitudes occur on a proper superset of the digraphs that allow restriction to the real numbers 

or, even further, the rationals. Consequently, among these fields, complex numbers guarantee the 
t***'* ' largest possible choice of topologies supporting a discrete quantum evolution. A similar construc- 

tion separates complex numbers from the skew field of quaternions. The result proposes a concrete 
ground for distinguishing between complex and quaternionic quantum mechanics. 

INTRODUCTION 

a 

Since the early 2000s substantial attention has been given to various types of quantum dynamics defined 
with respect to an underlying network, or graph. Studies have proposed such dynamics, mainly in the 
context of quantum information processing, to efficiently induce useful probability distributions. Relevant 
applications range from techniques for searching and exploring combinatorial structures 0, 0] to settings 
for universal computation [8j. Such models also play a central role in protocols for state preparation and 
transfer in nanodevices based on spin systems Q and in modeling the transport of energy in biochemical 
structures [24j . In most of these settings a graph specifies a symmetric network of interactions, but we need 
', to distinguish directions and will work instead with digraphs. 

' A digraph is denoted by G = (V, E) where V consists of n vertices and E is a set of ordered pairs of 

vertices, called directed edges. The source and target of an edge need not be distinct; we allow loops. 
(Equivalently, a digraph is a mathematical relation on a finite set.) A canonical way to associate a Hilbert 
space Hv to G is to define Hv = C™ = : i e V), where \i) is, as usual, the i-ih element of the standard 
basis. The space "Hv is the state space of a scalar quantum particle constrained to evolve on the vertices of 
G. To each digraph G is associated an adjacency matrix, which, over the class of digraphs we consider, may 
be any square matrix with entries in {0, 1}. For any vertex with a loop, the corresponding diagonal entry is 
1. To any matrix combinatorial matrix theory associates a zero-pattern [5], a matrix with entries in {0, ^} 
where ^ stands for an unknown entry that can take any value other than zero. (The term "sparsity pattern" 
is also used.) We adapt terminology introduced by Aaronson and Ambainis Q (see also Osborne [25|) to 
say that a unitary matrix U e U (n) (not necessarily Hermitian), or the quantum dynamics it specifies, is 
saturated Z-local on G if U has the same zero-pattern as the adjacency matrix of G. Supposing that H is 
a simple (non-directed) graph and F is the digraph on the same vertex set that has a symmetric pair of 
directed edges for each edge of H and that has a loop at each vertex, then a quantum dynamics on H is 
Z-local in the sense of Q if it is saturated Z-local on a subdigraph of F. 

The notion of Z-locality is applied to a discrete dynamics, defined for example by a quantum circuit model, 
just as the notion called i7-locality applies to evolution under a Hamiltonian. In fact, iJ-locality expresses 
locality by the interactions specified by the Hamiltonian operator (i.e., the clustering of correlations; see 
[25l]). Aspects of Z-locality for matrices in various stochastic ensembles have been investigated within the 
analysis of models describing the time evolution in quantum graphs and the quantum mechanics of systems 
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that are classically chaotic [l7, 2^, 3^]. Recall that a matrix M is unistochastic (resp. ortho stochastic) if 
there is a unitary (resp. real orthogonal) matrix U such that Mj ; j = |[w]i,j| 2 , for every i and j. Because of 
this definition, it appears obvious that Z-locality can be studied equivalently in the two settings of stochastic 
and quantum processes. 

Among other areas 12[ , these matrices occur in foundational questions [18 . [l9j and in high energy physics, 



where their role is to characterize interactions between elementary particles, with the Cabibbo-Kobayashi- 
Maskawa matrix [l6| as the arguably most famous example. 



Differences between stochastic ensembles have been discussed since the 1960s [23j and it is well-known 
that the set of unistochastic matrices includes properly the set of the orthostochastic ones. A complete 
characterization of saturated Z-locality is not currently known, although the matrix analysis literature 
contains several graph-theoretic and linear algebra conditions attempting to tackle this problem (see flpj ] 
and the list of references contained therein). Even if the original question was mathematically formulated in 
1988 the connection with minimum rank problems was underlined only recently 31 1. The connection is 



important because it locates the problem within the context of rank-related topological invariants including, 
for example, the Colin de Verdiere number and the Lovasz function [l3j]. 

While it is immediate to see that the zero-patterns of unitary (resp. real orthogonal) and unistochastic 
(resp. orthostochastic) matrices are identical, it is not obvious whether every unitary matrix can be "flat- 
tened" to an orthogonal matrix having the same zero-pattern. As a part of the effort to characterize the 
ortho- and unistochastic property [33], this question was resolved in the positive for matrices of dimension 



n ^ 5 [111 ], meaning that for small matrices the sets of orthogonal and of unitary zero-patterns coincide. It 
turns out that the general situation is different, as it is shown in the present work. 

For dynamics with amplitudes in a skew field F, we study in each case the set Z (F) of zero-patterns of 
digraphs that allow for saturated Z-local dynamics. We are interested specifically in the skew fields <Q>, M, C, 
and H of the rationals, the reals, the complex numbers, and the quaternions, respectively, whose sets of F- 
orthogonal matrices (preserving an abstract inner product over F) are, respectively, the rational orthogonal, 
orthogonal, unitary, and hyperunitary matrices. The set Z (F) consists of all zero-patterns of matrices that 
are F-orthogonal. We shall give a proof of the following statement: 

Theorem. Z (Q) c Z (R) c Z (C) c Z (H). 

A quantum dynamics on a graph G governed by a Z-local unitary matrix can be seen as a generalization 
of a (discrete) random walk on the vertices of G where the transition probabilities, which define a stochas- 



tic matrix, are substituted by transition amplitudes. Sinkhorn's theorem [21| states that every stochastic 
matrix can be made bistochastic without modifying its zero-pattern (recall that a nonnegative real matrix 
is stochastic if the entries of each row sum up to one; bistochastic if also the columns satisfy this property). 
However, we know that the zero- patterns of bistochastic and unistochastic matrices do not coincide [5|. This 
means that not every digraph G supporting a saturated random walk has a unitary matrix that is saturated 
Z-local on G. 

We remark that although we have defined Z (F) over the class of digraphs, the results of the theorem 
apply equally well to the more usual setting of symmetric networks. Each digraph G on n vertices, with 
adjacency matrix A, is associated uniquely to a bipartite simple graph B{G) on 2n vertices whose adjacency 
matrix has a block form with zero matrices in the diagonal blocks and A and its transpose in the off-diagonal 
blocks. The graph B{G) supports saturated Z-local dynamics over F if and only if G e Z (F). In addition, 
the dynamics of B(G) can without loss of generality be taken to be Hcrmitian. 

The theorem asserts that a characterization of potentially orthogonal digraphs must take the field into 
account. A corollary of the result is that there are digraphs, and indeed graphs, that require complex numbers 
to define a saturated Z-local quantum dynamics. In other words, complex numbers define saturated discrete 
quantum walks on a larger set of graphs than permitted by the reals; the same relationship holds between 



real and rational numbers. This last point solves the main open problem formulated in 11|. 

On the other hand, techniques have been proposed to define a quantum dynamics on every graph by the aid 
of some extra degrees of freedom that permit to enforce Z-locality on a topologically equivalent object. The 
most common techniques of this type are coined quantum walks, where the dynamics is lifted to a product 
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space composed by a shift and a coin register Q (or, equivalently, with the use of the graph-theoretic notions 



of a line digraph [27]); and Szegedy's generalization of Markov chains [29]. In this respect, the consequences 
of our theorem are significant when Z-locality is not reflected by the use of ad hoc constructions. More 
generally, Z-locality is not an obstacle from a point of view embraci ng spec ific algorithmic applications, 
because there are methods to translate between different number fields HI 22|. 

The next section gives a proof of our main result after stating a required lemma (proof in appendix). 
We describe in this work a concrete framework to distinguish between complex and quaternionic quantum 
mechanics. The debate about this topic is briefly addressed in the final section. 



PROOF OF THE THEOREM 



The three strict inclusions claimed in the theorem must arise in each case from a zero-pattern that belongs 
to one collection but not to the other, where the only allowed constraint in Z (F) is to make each entry 
zero or nonzero. Such a coarse degree of control makes it difficult to engineer an obvious dependence on the 
numerical field. We make use of a technique that allows us to impose an additional special type of constraint 
on a small zero-pattern with the guarantee that the special constraints are faithfully reflected by a pure 
zero-pattern that may be several times larger. 

Given a rectangular zero-pattern T, we define a four-way democracy on T to specify, within a single 
column of T, four ^-entries that are required to have the same magnitude. (The choice to define four-way 
democracies rather than another size is tied to a particular 7x7 zero-pattern whose rigidity properties are 
central to the proof of the lemma.) 

A pair (T, C), where T is a zero-pattern and C is a set of four- way democracies on T, is called a constrained 
zero-pattern. Given a constrained zero-pattern (T, C) and a skew field F, if a matrix A over F has zero- 
pattern T and has orthogonal columns we call A an orthogonal representation of T over F. If in addition 
the magnitudes of entries of A respect all the constraints of C, we say that A is a constrained orthogonal 
representation of {T,C) overF. 



Lemma. Let (T, C) be a constrained zero-pattern. Then there exists a zero-pattern T", containing 
T as a submatrix (so that (T',C) is also well-defined as a constrained zero-pattern), such that for any 
F e {Q, R, C, H} the following are satisfied: 

1. Every orthogonal representation of T' over F is also a constrained orthogonal representation of (T", C) 
over F. 



2. Every constrained orthogonal representation of (T, C) over F can be completed to an orthogonal rep- 
resentation of T" over F. 



3. If the constraints of C affect only a single column c of T, then column c in T' has no ^-entries other 
than the ^-entries of c in T. 



Letting (E, F) represent any consecutive pair in the sequence (Q, R, C, H), we exhibit a constrained zero- 
pattern (Tf,Cf) such that the zero-pattern T F promised by the lemma has an orthogonal representation over 
the skew field F but no orthogonal representation over the field E. To indicate which four ^-entries belong to 
a four-way democracy, we give them the same subscript (including multiple subscripts if a ^-entry belongs 
to more than one four- way democracy). With this notation, the three constrained zero-patterns (Tf,Cf) to 
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which we will apply the lemma are as follows: 
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First inclusion: Z (Q) cz Z(R). In the constrained zero-pattern (Tr,Cr) the five nonzero entries must all 
have the same magnitude. We apply the lemma to obtain a zero-pattern T^. Since the matrix A = [11111]^ 
is a constrained orthogonal representation of (Tr,Cr), it can be completed to an orthogonal representation 
A' of Tg. The orthogonal columns of A' can be scaled and completed to an orthonormal set of columns of a 
square matrix A", whose zero- pattern P belongs to Z (R). Now suppose by way of contradiction that P also 
belongs to Z (Q), implying the existence of a rational orthogonal matrix B", with zero-pattern P. Then a 
subset B' of the columns of B" is an orthogonal representation of Tj^ over Q whose columns are unit vectors, 
which by the lemma is also a constrained orthogonal representation of (T^,Cr) over Q. Since Cr applies to 
a single column, there are exactly five nonzero entries in that column of B', all with the same magnitude, 
namely l/v5, which cannot be the magnitude of a rational number. By contradiction, P does not belong to 
Z(Q), implying Z (Q) czZ(M). 

Second inclusion: Z (R) cz Z (C). The argument proceeds as in the proof of the first inclusion: (7c, Cc) 
has a constrained orthogonal representation consisting of two columns over C (whose inner product requires 
three complex numbers of the same magnitude to sum to zero) but does not have a constrained orthogonal 
representation over R. Thus T c has an orthogonal representation over C which can be scaled and completed 
to a unitary matrix of zero-pattern P e Z (C) , but P cannot belong to Z (R) because any orthogonal 
representation of T c over R would contain, as a submatrix, a constrained orthogonal representation of 7c 
over R. Thus we have Z (R) cz Z (C). 

Third inclusion: Z (C) cz Z(M). For the previous claim we used that fact that three numbers x, y, and 
z of equal magnitude cannot sum to zero in R, but can do so in C. In fact, if w\ and u>2 are the two 
non-real third roots of unity in C, we can have x + y + z = 0in precisely two ways: {y = uix, z = uJiy} or 
{y = u>2X,z = L02y}- Now we show that Th docs not have a constrained orthogonal representation over C. 
Suppose on the contrary that it did, and for the moment limit attention to just the first three rows of Th, 
naming the columns of that submatrix vo,vi,v 2 ,v 3 e C 3 . Since vi is orthogonal to v , three numbers of 
equal magnitude must sum to zero, with a choice to be made of u± or U2- If V2 is orthogonal to vo using 
the same choice, then v 2 will be parallel to vi; otherwise, since W2 = ujf, vi and v 2 will be orthogonal. In 
neither case can the remaining constraints be satisfied to construct a constrained orthogonal representation 
of the full pattern Th over C. Over H, however, there is a continuous family of third roots of unity, and Vi 
and V2 can be chosen so that both are orthogonal to Vo but so that they are neither parallel nor orthogonal 
to each other, and so that vi is orthogonal to v 3 but v 2 is not. This allows a complete construction of a 
constrained orthogonal representation of Te over H. It follows that some zero-pattern P, which completes 
the columns of Tj^ to a square matrix, separates the sets Z (C) cz Z (H). 



CONCLUSIONS 



Several questions of a mathematical nature remain open: What are the smallest examples of digraphs that 
distinguish between the number fields? What is the computational complexity for determining whether a 
given digraph admits a saturated Z-local unitary matrix? (In an email to the authors, S. Aaronson observed 
that the problem is reducible to the existential theory of reals and it is therefore in PSPACE.) Is there a 
combinatorial way to characterize these different families of digraphs? 

From the physical point of view, our result indicates a novel, alternative ground to explore the distinctions 
between real, complex (CQM), and quaternionic (QQM) quantum mechanics, a subject of debate whose 
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roots go back to the 1930s and the axiomatization of Birkhoff and Von Neumann (see [!,[28j]). In the study 
of correlations, quantum mechanics is "sandwiched" between classical mechanics and general probabilistic 
theories (see [7J for a recent treatment of this point in the graph theoretic framework). When we consider 
number fields, we seem to face a similar situation: while the choice of number field does not affect the 
computational power of the theory [TR [22j , there are cogent arguments about the inadequacy of real numbers 
{e.g., parameter counting for bipartite mixed states, continuity of time, the quantum de Finetti theorem, 
the need of superselection rules, etc. [l[). 

In the other direction, the status of the connection between CQM and QQM is still unresolved. The 
combinatorial tool of saturated Z-local dynamics may contribute to the separation of these theories in terms 
of multiparty correlations obtainable with given resources (such as a fixed Hilbert space dimensionality 
specified according to the physical system under consideration). It is possible at the experimental level that 
observable entanglement measures could be employed to reach a contradiction (in terms coherent with the 
known theory) to show that QQM leads to provably unphysical claims or that it has a particular range of 
applicability. The spirit of these latter assertions is clearly speculative and they suggest directions where 
further work is required. 
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APPENDIX: PROOF OF LEMMA 



The lemma used in the proof of our main result asserts that one can translate the constraints of a collection 
of four-way democracies on a small zero-pattern matrix into the zero-pattern alone of a larger matrix. The 
mapping of constraints depends on a special property of the following 7x7 matrix, 



M 



1 

1 -1 



1 -1 
1 1 

1 

1 
1 



and its zero-pattern, which we call F. It will be convenient to index the rows and columns of 7 x 7 matrices 
by the sequence (0, 1, ... , 6) and to stipulate that index arithmetic is modulo 7. We first observe a few facts 
that are interesting but not directly essential to the proof of the lemma: 

1. Entry (i, j) of F is * if and only if i — j is a quadratic residue modulo 7. 

2. The zero-entries of F give incidence relations on the points and lines of a Fano plane. 

3. The ^-entries of F give the negated entries in a multiplication table of the imaginary part of the 
octonions. 

4. M + I is the adjacency matrix of a tournament graph. 

5. The matrix is rational orthogonal. 



G 



The last two points are interesting partly because it is an open problem [20J whether there exists an orthogonal 
(or indeed unitary) matrix of size greater than 3x3 that is saturated Z-local on a tournament graph. The 
heuristic of dimension counting is pessimistic about the existence of such a matrix: Since 0(n, M) is a 
manifold of dimension n(n — l)/2 and each zero-entry introduces one constraint, non-redundant constraints 
would give, for the set of orthogonal matrices that are saturated Z-local on any given tournament graph, 
a supposed dimension that is less (by n) than zero-dimensional, or in other words a decidedly empty set. 
Observe, however, that if a loop is added to each vertex of a tournament graph the expected dimension is 
now zero, so that one would generically expect a finite set of rigid solutions. Such a looped tournament 
graph is the form that F takes, and the special property we will establish for M and F is indeed a form of 
rigidity. 

Claim. Suppose that A is an orthogonal representation of F over F, where F e {Q,K, C,H}. Then there 
exist 7x7 diagonal matrices U and D over F, with U additionally hyperunitary, such that UAD = M . In 
particular, A is a constrained orthogonal representation of (F, C) over F, where C consists of one four-way 
democracy on the ^-entries of each column. 

Proof of claim. We construct the matrices D and U in stages, D as a product of seven matrices D = 
DqDi . . . De and U as a product of six matrices U = U§ . . . U2U1, where each Di or Ui may differ from 
the identity matrix in the entry but nowhere else. We also recursively define a sequence of matrices 
j4o, . . . , Aq and Bq, . . . , B$ in a way that uniquely specifies, for any A, the choice of each Di or Uf. 

• We start with A = A. 

• For i = 0, . . . , 6, we choose Di such that Bi = AiDi has entry (i, i) equal to —1. 

• For i = 1, . . . , 6, we choose Ui such that Ai = f/j_Bj_i has entry (i, i — 1) real and positive. 

The final result of this is B = Bq = U AD, which (like M) has —1 along the diagonal, and which has positive 
real entries on the subdiagonal. (Note that the entire process can be completed without leaving the field F.) 
We name the entries of B as follows: 
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-1 
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x 5 
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UAD 



The entries xq through x§ are by construction real and positive but xq might, a priori, have a non-trivial com- 
plex or quaternionic phase. Since U is unitary and both A^A and D are diagonal, the product D* A'U'UAD 
is also diagonal and the columns of B are mutually orthogonal. By cyclic symmetry, any equation derived 
from the orthogonality of columns of B must remain true when the set of indices is permuted cyclically. We 
introduce the shorthand X$ = (The conjugation in this definition is only necessary in the case of xq, 

whose quaternionic phase is for the moment undetermined.) The fourteen variables {xi,Xi} all commute, 
even in the case F = H, since at most one related pair is not real. 

From the fact that entry (1,0) of B^B is zero, we conclude that xiyo = xq and thus, for all i, 

yi = XiXi+i. 

From the fact that entry (2,0) of B^B is zero, we conclude that y^ZQ = yo and thus, for all i, 

Zi = XiXi + iXi + 2Xi + 3- 

From the fact that entry (4,0) of B*B is zero, we conclude that = zq and thus, for all i, 

Xi = XiXi + iXi + 2Xi + 3Xi + 4Xi + 5Xi + gXi + j 
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and, for example, 



X XiX2X 3 X4X 5 X 6 = 1, 

which in particular shows that x 6 is after all real and positive, so that in all cases = x" 1 . It is convenient 
to take logarithms of this equation, together with its cyclic index permutations, and write the homogeneous 
system of linear equations 
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which has unique solution {lnx^ = 0}. This implies UTD = M and completes the proof of the claim. 

We have in F an example of a zero-pattern which is unitarily rigid, meaning that there is a unique unitary 
matrix with pattern F up to scaling of rows and columns by complex phase. A matrix analogous to \M can 
be constructed of size p x p for other primes p = 4fc — 1 (for k an integer), putting (1 — k)/k on the diagonal 
and 1/k in the off-diagonal positions where the ordered index difference is a quadratic residue modulo p, 
with zeros elsewhere. In the cases k = 1 and k = 2 the resulting zero-pattern is unitarily rigid, and it is a 
natural question, to which we do not know the answer, whether this gives a unitarily rigid zero-pattern for 
every prime congruent to 3 modulo 4. 

The rigidity of F is what allows us to construct large zero-patterns from small constrained zero-patterns. 
Consider the case of the constrained zero-pattern (T, C) where T is the 7x7 pattern F and where C consists 
of four-way democracies on any subset of the columns. In this special case, T" = T already satisfies the 
conclusions of the lemma. In general we will construct a zero-pattern T" that positions several copies of F 
as submatrices (one for each element of C), in each case with one column overlapping the submatrix T and 
the other six columns outside of it. We give the name F 6 to the 7x6 zero-pattern consisting of the last six 
columns of F. 

Proof of lemma. Let (T, C) be a constrained zero-pattern. It suffices to prove a special case of the lemma in 
which C is a singleton {Co}, since the process can be iterated: If C = {Co, . . . , C n _i} then we first apply the 
special case of the lemma to (T, {Co}) to obtain T{, then apply the special case of the lemma to (T[, {Ci}), 
to obtain T 2 , and so forth until we obtain T' n = T". Assuming then that C = {Co}, we start with T and 
append rows or columns in three stages to obtain R, S, and finally T' . 

First, we obtain the zero-pattern R from T by appending three rows of all-zeros. Call these the control 
rows. 

Next, we obtain the zero-pattern S from R by appending six new columns, the control columns, whose 
zero-pattern will depend on the six columns of the pattern F$: The four rows that correspond to Co get the 
pattern of rows 0, 1, 2, and 4 of Fq, and the three control rows get the zero-pattern of the remaining rows 
3, 5, and 6 of Fq. All other rows are zero. Now we can say some things about orthogonal representations 
of S. Let the column of T or of S that is constrained by Co be called Co, and call the control columns 
ci, . . . , eg. Other than rows containing at most a single the zero-pattern of Co and the control columns is 
F. It follows that any orthogonal representation of S over F is also a constrained orthogonal representation 
of (S, {Co}) over F, so that S already fulfills the first conclusion of the lemma. However, there may still 
be constrained orthogonal representations of (T, {Co}) over F that cannot be completed to an orthogonal 
representation of 5* over F. This can happen because the ^-entries of some non-co column a of T and the 
^-entries of some control column b may overlap in as many as two rows. The purpose of the third stage is 
to remedy any resulting non-orthogonality. 

Finally, we obtain the zero-pattern X" from S by appending several additional rows, in each of which 
exactly two entries are * and the rest are zero. For every pair (a, b) consisting of a column a ^ Co from the 
original zero-pattern T and a column b from the control columns, we append either zero, one, or two rows 
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with ^-entries at a and b. The mutual support of a and b is the set of rows in S where column a and column 
b both take the value ^; this can either be empty or can consist of one or two rows. If the mutual support 
of a and b consists of exactly one row, then for the pair (a, b) we append exactly one new row while building 
T' . The two new ^-entries in this row have no other constraints and only ever interact with each other, and 
so columns a and b, which in S were forced to be non-orthogonal, can now easily be made orthogonal in T'. 

If the mutual support of a and b is empty, then columns o and b are automatically orthogonal, and we 
append no additional row for the pair (a, b). 

If the mutual support of a and b consists of two rows, then regardless of S we may, to ensure that columns 
a and b can easily be made orthogonal, always safely append two rows for the pair (a, b) when building T' . 
(If there is reason to prefer a smaller T", in most cases inspection is sufficient to assure that appending no 
additional row for the pair (a, b) will serve just as well.) 

Each appended row or pair of rows in T' allows for the orthogonality of a pair of columns a and b 
without spoiling the orthogonality previously secured for any other pair of columns. The complete process 
of appending rows ensures that every constrained orthogonal representation of (T, {Co}) over F can be 
completed to an orthogonal representation of T' over F. The conclusion also remains, as it did for S, that 
the zero pattern F occuring as a submatrix forces every orthogonal representation of T" over F to be a 
constrained orthogonal representation of (T", {Co}) over F. Since column c is always distinct from a and 6, 
it acquires no new ^-entries in this process, which means, as the process is iterated, that if all the constraints 
in C affect only a single column c, then column c in T" has no ^-entries other than the ^-entries of c in T. 
This completes the proof of the lemma. 
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